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Abstract In the current paradigm for the modulation of galactic cosmic rays (GCRs), dif- 
fusion is taken to be the dominant process during solar maxima while drift dominates at 
minima. Observations during the recent solar minimum challenge the pre-eminence of drift 
at such times. In 2009, the ~2 GV GCR intensity measured by the Newark neutron monitor 
increased by ~5% relative to its maximum value two cycles earlier even though the average 
tilt angle in 2009 was slightly larger than that in 1986 (~20° vs. ~14°), while solar wind 
B was significantly lower (~3.9 nT vs. ~5.4 nT). A decomposition of the solar wind into 
high-speed streams, slow solar wind, and coronal mass ejections (CMEs; including post- 
shock flows) reveals that the Sun transmits its message of changing magnetic field (diffusion 
coefficient) to the heliosphere primarily through CMEs at solar maximum and high-speed 
streams at solar minimum. Long-term reconstructions of solar wind B are in general agree- 
ment for the ~ 1900-present interval and can be used to reliably estimate GCR intensity over 
this period. For earlier epochs, however, a recent l0 Be-bascd reconstruction covering the 
past ~10 4 years shows nine abrupt and relatively short-lived drops of B to <0 nT, with the 
first of these corresponding to the Sporer minimum. Such dips arc at variance with the recent 
suggestion that B has a minimum or floor value of ~2.8 nT. A floor in solar wind B implies 
a ceiling in the GCR intensity (a permanent modulation of the local interstellar spectrum) at 
a given energy/rigidity. The 30-40% increase in the intensity of 2.5 GV electrons observed 
by Ulysses during the recent solar minimum raises an interesting paradox that will need to 
be resolved. 
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1 Introduction 

Forbush (1954, 1958) discovered the 11 -yr galactic -cosmic-ray (GCR) modulation cycle 
and noted its anti-correlation with the sunspot cycle. Parker (1965), Gleeson and Axford 
(1967), Jokipii and Parker (1970), and others [see Moraal 1976, 201 1 (this issue), Potgieter 
1998, and Heber and Potgieter 2006 for reviews] developed the particle transport equation 
which has terms governing GCR diffusion, convection, curvature/gradient drift, and adia- 
batic deceleration in the heliosphere. During the past three decades, modulation research has 
been primarily focused on understanding and specifying the roles of two of these terms — 
diffusion and curvature/gradient drift: 

(1) In the diffusion-centered view, modulation results from the outward propagation, at the 
solar wind speed, of large steplike decreases in the GCR intensity (McDonald et al. 
1981; cf., Morrison 1956, and Lockwood 1960). The disturbances responsible for these 
decreases came to be called “merged interaction regions” (MIRs) or “global merged 
interaction regions” (GMIRs; Burlaga et al. 1984, 1985, 1993; McDonald ct al. 1993). 

(2) Alternatively, the approach championed by Jokipii and colleagues (Jokipii et al. 1977; 
Jokipii and Thomas 198 1 ; Kota and Jokipii 1983) stresses the importance of curvature 
and gradient drifts of GCRs in the large-scale field of the heliosphere. In an early success 
of the drift picture, Jokipii and Thomas (1981 ) predicted the peaked cosmic ray intensity 
maximum in 1987 during the solar minimum between cycles 21 and 22 (Smith 1990). 

Today, these two approaches — one emphasizing propagating diffusive barriers (e.g., Perko 
and Fisk 1983; Burlaga et al. 1985) and the other focusing on particle drift effects as para- 
meterized by the solar “tilt” angle (e.g., Jokipii and Thomas 1981 ; Smith and Thomas 1986; 
Le Roux and Potgieter 1992a, 1992b) — have been combined to form the current paradigm 
for GCR modulation (e.g., Jokipii and Wibberenz 1998; Potgieter 2008, 2010). In this pic- 
ture, drift effects are dominant on the decline of the solar cycle and at solar minimum while 
diffusion is dominant at solar maximum. 

The current picture of modulation at solar maximum was challenged by Cane et al. 
(1999) and Wibberenz et al. (2002) who argued that since steplike decreases in the GCR 
intensity were first observed at Earth and then later in the outer heliosphere, it was incon- 
gruous to attribute them to MIRs or GMIRs formed by the merger/coalescence of coronal 
mass ejections (CMEs), shocks, and high-speed streams between 10-15 AU. Instead, they 
attributed modulation to global changes in the Sun’s magnetic field. Subsequently, Sheeley 
et al. (2001), Cliver and Ling (2001b), and Wang et al. (2006) focused on changes of the 
Sun’s low-latitude or non-axisymmetric open flux, which are more closely related to the 
sunspot number and solar wind B, than is the total open flux. More recently, the behavior 
of the Sun and the related cosmic ray response during the extended solar minimum period 
following cycle 23 have raised questions about the dominance of drift at this phase of the 
cycle. 

In this study, we reinvestigate the solar and solar wind drivers of GCR modulation, both 
for 1 1-yr and longer-term modulation. Our focus in this paper is on the modulation of high- 
rigidity (~2 GV) positively charged particles at 1 AU. Our analysis is presented in Sect. 2 
and the results are summarized and discussed in Sect. 3. 
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2 Analysis 

2. 1 Diffusion vs. Drift 

Figure ! contains timelines of solar wind B , the cosmic ray intensity (inverted) measured by 
the Newark neutron monitor (cutoff rigidity of ~2 GV), and the sunspot number for solar 
cycles 20-23. In drift formalism, the intervals between maxima of the cycles are designated 
as A positive or A negative where A positive (negative) cycles correspond to those for which 
the Sun’s dipole magnetic field points outward (inward) in the northern hemisphere. During 
epochs when A > 0, positively charged GCRs approach the inner heliosphere from over 
the poles (Jokipii and Kopriva 1979). In epochs when A < 0, the preferred direction of 
approach of positively charged GCRs is inward along the heliospheric current sheet (HCS) 
that separates oppositely-directed fields in the northern and southern hemispheres of the 
heliosphere. The average excursion of the HCS from the heliographic equator is referred to 
as the tilt angle (a) of the heliospheric current sheet (available at http://wso.stanford.edu/). 
In practice a is defined as the average of the maximum angular excursions of the coronal 
neutral line (the base of the HCS) north and south of the Sun’s equator during a 27-day solar 
rotation. The tilt angle is widely used to characterize drift effects on cosmic ray intensity, 
and we shall do so in this study. It should be kept in mind, however, that tilt angles provide 
only a crude approximation of the magnetic architecture of the heliosphere and that more 
sophisticated parameters, such as the total length of the wavy current sheet between Earth 
and the outer modulation boundary (e.g., Le Roux and Potgieter 1990) are needed. 

Figure 2 contains scatter plots of 27-day averages of the GCR intensity measured at 
Newark vs. the tilt angle for A positive (a) and A negative (b) epochs, from 1976 to 2009, 
the start of cycle 21 to the end of cycle 23. In each plot, blue filled points indicate data 



Fig. 1 Solar wind B, Newark cosmic ray intensity (inverted), and sunspot number, 1964-2009. The Sun’s A 
positive and A negative magnetic epochs are indicated 
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Fig. 2 Scatter plots, with regression lines, of 27-day averages of the Newark neutron monitor counting rate 
vs. the “classic” tilt angle (which uses a line-of-sight boundary condition at the photosphere and includes a 
significant polar field correction) of the HCS from Wilcox Solar Observatory for (a) A positive and (b) A 
negative polarity epochs. Individual rise and fall intervals of the solar cycle are color coded. Intervals when 
a > 70° (undefined tilt angle) have been excluded from the plots 


from cycle 21, black filled points from cycle 22, red filled points from cycle 23, and blue 
open points from cycle 24. Thus, the rising phase of cycle 21, when A > 0, is shown by the 
solid blue symbols in (a), while data for the declining phase, when A < 0, are shown in (b). 
Note that the GCR intensity is a multi-valued function of the tilt angle for both polarity 
epochs, i.e., a single tilt angle corresponds to more than one GCR intensity. For example, 
during the A positive epoch in the 1990s, the recovery phase of modulation (black points 
in (a)) exhibited a much stronger dependence on tilt angle than did the onset phase (red 
points). The multi-valuedness of the functions in both Fig. 2(a) and (b) partially reflects the 
inadequacy of the tilt angle as a parameter for the drift effect. On the rise of a solar cycle, the 
large tilt angles in the plots (~60-70°) will correspond to a shorter total length of the HCS 
than is the case for the decay phase. Thus, all else being equal, one would expect modulation 
to be greater for high tilt angles during the recovery of GCR intensity (resulting in the black 
points lying below the red points in (a)). Note in (b), however, that for low till angles, the red 
points corresponding to the decay of cycle 23 lie mainly above the solid blue points for the 
decay of cycle 21 while the blue circle data points for the rise of cycle 24 lie above the black 
points for the rise of cycle 22. In these cases, the multi-valuedness of the GCR intensity vs. 
tilt relationship suggests that some factor other than drift is driving the modulation. We will 
pursue this point in greater detail below. 

Figure 3 shows the Newark GCR intensity plotted against solar wind B (27-day aver- 
ages). In this case, the GCR intensity is more nearly a single-valued function, albeit with 
significant scatter, of solar wind B for the rise and fall phases within a given polarity epoch. 
However, there is a clear difference in the response of the GCR intensity to changes in B in 
A positive and negative epochs. For low B values (< 6 nT) in A positive epochs (Fig. 3(a)), 
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Fig. 3 Scatter plots, with regression lines, of 27 -day averages of the Newark neutron monitor counting rate 
vs. 27-day averages of solar wind B for (a) A positive and (b) A negative polarity epochs. Individual rise and 
fall intervals of the solar cycle are color coded. Rise portions include a > 70° intervals. The regression line 
in (a) is a second order fit 

GCR intensity decreases more slowly as B increases (cf., Wibberenz et al. 2002) than is the 
case for A negative epochs (Fig. 3(b)). 

This effect can also be seen in Figs. 4 and 5. In Fig. 4(a) and (b), we consider the GCR 
response to changes in B for fixed ranges of the tilt angle in negative and positive polarity 
epochs, respectively. Note that in A negative epochs (Fig. 4(a)), the slope is comparable 
over the full range of tilt angles. In contrast, during A positive epochs (Fig. 4(b)), the re- 
sponse for low tilt angles is weak (neglecting the circled outlier in the correlation analysis) 
but increases rapidly with increasing tilt angle. Figure 4 is a direct demonstration of the 
suggestion by Wibberenz et al. (2001) that “for A > 0 polarity epochs, the sensitivity of the 
CR response to variations in B is partly controlled by the size of the tilt angle, at.” Figure 5 
contains time lines of: (a) the tilt angle and Newark cosmic ray intensity (inverted, 27-day 
averages) for 1976-2009, and (b) solar wind B (27-day averages) and Newark GCR inten- 
sity (inverted, 27-day averages) from 1965-2009. While the correlation coefficients between 
the parameters in the two panels are not greatly different (e.g., r =0,80 in (a) vs. 0.74 in 
(b) for 27-day averages from 1976-2009), the closer relationship of B with the cosmic ray 
intensity for the common years in the plots (in comparison with the relationship between a 
and the cosmic ray intensity), both for phase and amplitude, is apparent. The weaker fidelity 
between the tilt angle and GCR intensity is noticeable on the rise of the odd-numbered solar 
cycles 21 and 23, an effect reported previously by Cliver and Ling (2001a) (see also Na- 
gashima and Morishita 1980, and van Allen 2000), and is particularly evident during the 
decline of cycle 23 and the rise of cycle 24. It is instructive to bring the data in Fig. 5 to 
bear on the multi-valued nature of the relationship between the GCR intensity and the tilt 
angle in Fig. 2. Specifically, we note: (i) the large separation between the red and black 
points for tilt angles of ~60° in (a) corresponds to a ~2.5 nT difference in solar wind B 
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Fig. 4 (a) Scatter plots of the Newark neutron monitor counting rate (or cosmic ray intensity (CRI)} vs. solar 
wind B for three ranges of the tilt angle during A negative polarity epochs. The color coding is the same as 
in Figs. 3 and 4. (b) Scatter plots of the Newark CRI vs. solar wind B for three ranges of the tilt angle during 
,4 positive polarity epochs. The circled outlier (Carrington Rotation 1856; 20 May-16 June 1992) in the top 
panel was not used in the correlation analysis 


between 1998.5-1999.5 and 1991; (ii) the separation between the red and solid blue points 
for tilt angles of ~30° in (b) corresponds to a ~2.5 nT difference in solar wind B between 
2007 and 1984.5-1985.5; (3) and the separation between the blue circle and black points 
for tilt angles of ~45° in (b) corresponds to a ~2.5 nT difference in solar wind B between 
~2010 and 1988. In each case, the year with the larger value of B had the lower associated 
GCR intensity, consistent with the notion that solar wind B is the primary driver of GCR 
modulation at all phases of the 1 1-yr cycle. 

Figure 6 provides support for the idea that tilt angle variation does not drive modulation 
at solar minimum as is commonly thought to be the case (e.g., Potgieter and Le Roux 1994; 
Jokipii and Wibberenz 1998; Ferreira and Potgieter 2004). Figure 6(a) shows that while the 
lowest value of the annually-averaged tilt angle was essentially constant (~ 15-20°) at the 
last four solar minima (actually increasing slightly in 2009 relative to the three previous 
minima), the corresponding cosmic ray intensity was significantly (~5%) higher during the 
minimum in 2009. As shown in Fig. 6(b), the marked increase in cosmic ray intensity at 
Newark during this minimum was accompanied by a significant decrease (from ~5.4 nT to 
~3.9 nT) in solar wind B, relative to the three previous minima. 

A final point on the role of drift on the 1 1-year variation of GCR intensity — in ba- 
sic drift theory (Jokipii and Thomas 1981; Kota and Jokipii 1983) the tilt angle has the 
same or at least a non- A dependent evolution from cycle-to-cycle. Figure 7 shows that 
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Fig. 5 (a) Tilt angle of the HCS {red) and Newark neutron monitor counting rate ( black ; 27-day averaged, 
inverted, arbitrary units) from 1976-2009. (b) Solar wind B (red; 27-day averaged) and Newark neutron 
monitor counting rate (black; 27-day averaged; inverted, arbitrary units) from 1976-2009 


while that this assumption may be true for the rise of the cycle, it may not be true 
for the decline (Cliver 1993; Cliver et al. 1996; Cliver and Ling 2001a), In cycle 22, 
the tilt angle declined more rapidly in comparison with cycles 21 and 23. A relatively 
rapid decrease also appears to have occurred during cycle 20 (Hundhausen et al. 1980; 
Cliver 1993). A longer run of observations is required to establish the reality /regularity of 
this cycle-to-cycle variation in the tilt angle and to assess any associated systematic effects 
on the shapes of modulation cycles. The tilt angle record through 2010 indicates that thus far 
cycle 24 is evolving differently from the previous three cycles, with tilt angles ~45 0 reached 
more rapidly than in the past. 

2.2 How Does the Sun Transmit Its Message of Increased Magnetic Field to the 
Heliosphere? 

To answer this question, we will make use of a decomposition of the solar wind into its three 
principal components — slow solar wind (SSW), high-speed streams (HSSs), and transient 
structures including CMEs and post-shock flows, and an additional category for intervals 
when the flow type could not be determined, e.g., due to data gaps. This decomposition was 
made by Richardson et al. (2000, 2001, 2002) for the years 1972-2000 and has been updated 
through 2009, and extended back to 1965, for the present study. The decomposition is based 
on OMNI solar wind data and ancillary data from solar energetic particle detectors and 
neutron monitors, as discussed in Richardson et al. (2000). Figure 8 shows the contribution 
of each of the three flow types to annual averages of solar wind B for 1965-2009, with 
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Fig. 6 (a) Yearly-averages of Newark neutron monitor counting rate vs. yearly-averages of the tilt angle of 
the heliospheric for the last four solar minima, (b ) Yearly-averages of Newark neutron monitor counting rate 
vs. yearly-average of solar wind B for the last four solar minima. The minima are years with lowest average 
B and a values; for the 1986 and 1996 minima, the cosmic ray maxima occurred the following year, being 
higher by 15ctss _1 and 17 cts s^ 1 , respectively, than the plotted values 



Years from “Onset" of Solar Cycle 


Fig. 7 Superposed epoch plot showing the evolution of the tilt angle of the HCS during each of the last three 
solar cycles (updated from Cliver and Ling 2001a). The zero epoch was adjusted to obtain best agreement for 
the rise phases of the three cycles 
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Year 

Fig. 8 Cumulative distribution graph showing the contribution of SSW, HSSs, and CMEs to average B from 
1965-2009 (adapted from Cliver and Ling 2010). The dashed line is drawn at 2.8 nT, the revised position of 
the “floor” in the solar wind 


the contribution from the unclassified intervals distributed proportionately (according to the 
time Earth was immersed in SSW, HSSs, and CMEs during the year) among the three basic 
types. Note that the y-axis in Fig. 8 is not the average B in each flow type, although the sum 
of the three contributions adds up to average B for each year. Rather it is what average B 
would have been for a given flow type in the absence of the other two types (setting their 
contributions equal to zero and dividing the sum of hourly average B values in, e.g., CME 
flows, by the total number of hours of observations during a given year). We will refer to 
these contributions to B as Z?ssw> #hss> and #cme- 

Figure 8 shows that B ssw does not vary strongly with the solar cycle; it has a mean value 
of ~2.2 ± 0.4 nT over the 44-yr interval, B m $ peaks on the decline of the cycle in accord 
with the well-known tendencies for coronal holes and recurrent geomagnetic storms to be 
prominent during that phase of the cycle (e.g., Crooker and Cliver 1994). It appears that B 
from HSSs dominated modulation events (“minicycles”) on the decline of the solar cycle 
in 1974 (Hewish 2001; Cliver et al. 2003), 1984, and 1994, when Bq me is relatively small. 
Note that at the solar minima in 1976, 1986, 1996, and 2009, CMEs make a relatively small 
contribution to B. Thus the ~5 % increase in the Newark neutron monitor counting rate in 
2009 above the average for the three previous minima (Fig. 6) is attributed to the decrease 
in B hss during this year. The decrease in t9 H ss from ~2.8 nT during 1976/1986/1996 to 
~1.2 nT in 2009 accounts for the ~1.5 nT drop in B during these periods (from ~5.4 nT 
to ~3.9 nT). The contribution to B from CMEs tracks the solar cycle (Webb and Howard 
1994), and essentially disappears at solar minimum. Figure 8 shows that CMEs are primarily 
responsible for increasing B at solar maximum, and that HSSs are principally responsible 
for B variation at solar minimum as well as during modulation minicycles on the approach 
to minimum. Figure 9(a, b, c) shows scatter plots of GCR intensity vs. Z/ C me, B tiSS , and 
B ssw, respectively. The strongest anti-correlation is for CMEs; in fact, for SSW, there is a 
weak positive correlation. 
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Fig. 9 Scatter plots, with regression lines, of New ark GCR intensity vs. (a) Bq me> (b) 5 hSS< and (c) B ssw 

2.3 A Paradox Involving GCR Electrons at the Recent Solar Minimum 

Heber et al. (2009) noted that the modulation amplitude of GCR electrons during 2008 
exceeded the level reached during the mid-1990s minimum by 30-40%. In comparison, 
the modulation amplitude for the Newark neutron monitor in 2009 increased by only ~4% 
over the 1997 minimum. Modulation models predict that at true minimum the normalized 
count rates of both GCR electrons and protons should be equal (Potgieter et al. 2001). The 
failure of the modulation amplitude of the protons to reach a value commensurate with that 
for electrons is attributed (Heber 2010) to the fact that the tilt angle dropped below 15° 
for only 3 consecutive rotations in 2009, vs. 8 in 1995-1996. The lowest yearly (for any 
13 consecutive rotations) average of the tilt angle during the mid-1990s was ~13° [August 
1995-June 1996] vs. 19° for the corresponding interval during the first decade of the 21st 
century [December 2008-November 2009], Correlations based on cycles 21-23 (Fig. 10) 
show that the decrease of ~1.2 nT between the average values of B for 1997 and 2009 is 
far too small to account for a 30-40% increase in the modulation amplitude of GCR protons 
at Newark (a) and that the corresponding increase of ~5° in the tilt angle is far too small to 
negate such an increase (b). For example, as shown in Fig. 10(b), the average response of 
the Newark neutron monitor to a 10° change in the tilt angle is only ~3%. 

2.4 Long-Term Modulation 

2.4. 1 Emerging Consensus on Long-Term Reconstructions of B 

Long-term reconstructions of the solar wind magnetic field, once contentious (e.g., Lock- 
wood et al. 2006; Svalgaard and Cliver 2006; McCracken 2007), have recently reached a 
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Fig. 10 (a) Scatter plot (with regression line) of yearly values of the GCR modulation amplitude (MA) 
measured by the Newark neutron monitor from 1976-2009 vs. solar wind B (1987 = 100%). (b) Scatter plot 
(with regression line) of yearly values of the GCR modulation amplitude measured by the Newark neutron 
monitor from 1976-2009 vs. the tilt angle of the HCS ( 1987 = 100%) 


consensus for the 20th century. In Fig. 12 of their recent paper, Svalgaard and Cliver (2010) 
show the evolution and, for the most part, resolution between the discordant geomagnetic- 
based reconstructions of Lockwood et al. (1999) and Svalgaard and Cliver (2005). Further- 
more, in their Fig. 14, they show that the Steinhilber et al. (2010; cf., Webber and Higbie 
2010a, 2010b; Webber et al. 2010) 10 Be-based reconstruction of B, which superseded and 
extended that of McCracken (2007), is now in general agreement with the geomagnetic- 
based reconstructions for the last ~100 years. 

In Fig. 11, we use the regression relationship between the modulation amplitude of 
the GCRs measured by the Newark neutron monitor and annual averages of solar wind B 
(r = 0.89) from Fig. 10(a) and the reconstruction of B from 1835-2009 obtained by Sval- 
gaard and Cliver (2010) to reconstruct the modulation amplitude of ~2 GV protons for this 
interval and compare it with observations from 1963-present. Despite the gross simplicity 
of this approach, the reconstructed GCR modulation amplitude matches observations fairly 
well, except for cycle 20 (indicated by the blue oval). The good agreement within the oval 
between the reconstructions of the amplitudes based on actual B and on IDV-bascd B (from 
Svalgaard and Cliver 2010) indicates that the measurements of B during the early space age 
are sound and that the cause of the cycle 20 anomaly does not lie with erroneous B -values 
(as suggested, e.g., Belov 2000). Wibberenz et al. (2002) were able to model cycle 20 with 
reasonable fidelity using the observed B and the longer GCR recovery time (T ~ 100 days) 
characteristic of A negative cycles and solar maximum intervals. 
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Fig. 11 Reconstruction of annual values of the GCR modulation amplitude measured by the Newark neutron 
monitor based on the regression relationship in Fig. 1 0(a) and derived using: (a) observed solar wind B (blue 
line), and (b) IDV09-based B ( red line) from Svalgaard and Cliver ( 2010 ). The observed ~2 GV modulation 
amplitude ( black line) is shown for comparison. The dashed blue oval highlights the anomalous behavior of 
solar cycle 20 and the dashed black line indicates the ceiling in the cosmic ray intensity corresponding to a 
floor of ~2,8 nT 


2.4.2 A Floor in the Solar Wind Magnetic Field: Implications for GCR Modulation 

Recently, Cliver and Ling (2010) have argued that a floor exists in the solar wind magnetic 
field strength of ~2.8 nT (cf., Svalgaard and Cliver 2007). They attribute the floor to a con- 
stant baseline open flux of ~8 x 10 B Wb associated with the slow solar wind. Cliver and 
Ling suggest that the floor corresponds to a ground-state of the solar wind, achieved in the 
absence, or near absence, of the solar cycle, i.e., during periods such as the Maunder Mini- 
mum (Eddy 1976). Figure 12 shows that the existence of a floor in B implies a floor in the 
modulation potential ( <t») of MOO MV. Moreover, the regression equation from Fig. 10(a) 
and dashed line in Fig. 1 1 show that a floor in B of ~2.8 nT (or <J> of MOO MV) implies a 
ceiling in the cosmic ray intensity (M08% for Newark; 1987 = 100%). In a sense, the floor 
results in a permanent modulation of the local interstellar spectrum at any given energy or 
rigidity. 

The reconstruction of solar wind B by Steinhilber et al. (2010) extends back over the last 
MO 3 4 years (Fig. 13). During this interval there are nine periods, beginning with the Sporer 
Minimum (Eddy 1976) from ~ 1420—1570, where B dropped to zero, violating the notion 
of a floor. We will discuss these extreme dips further in Sect. 3. 


3 Conclusion 

3.1 Summary 

In this study of solar and solar wind drivers of 1 1-yr and long-term modulation of the ~2 GV 
galactic cosmic ray proton intensity, we have found the following: 
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Fig. 12 Solar wind B as a function of the modulation potential (after Steinhilber et al. 2010); the 
4> values were taken from: http://www.faa.gov/data_research/research/med_huinanfacs/aeromedical/media/ 
cari6mex.zip (see Fig. 2 in Steinhilber et al. 2010). A floor of in B of ~2.8 nT implies a minimum <t> of 
~ 100 MV 
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Fig. 13 Reconstruction of B (40-yr running means) based on 1() Be concentration in ice cores for the last 
~10 4 years (adapted from Steinhilber et al. 2010). The dashed green line drawn at ~4.6 nT represents the 
value of the floor reported by Svalgaard and Cliver (2007) and the dashed red line drawn at ~2.8 nT indicates 
the revised floor (Cliver and Ling 2010). The red ovals encompass sharp drops to B-values <0 nT 


(1) GCR intensity is a multi-valued function of the tilt angle of the heliospheric current 
sheet, both for different phases (rise and fall) of the solar cycle within a given A posi- 
tive epoch, and for the same phase of different solar cycles within A negative epochs. 
This multi-valuedness is due, in part, to the inadequacy of the tilt angle as a parameter 
for drift modulation. 

(2) GCR intensity is roughly a single-valued function of solar wind B within a given po- 
larity epoch, but the function differs between A positive and A negative epochs. In 
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particular, GCR intensity has a weaker response to changes in B near the onset of the 
modulation cycle (i.e., for low B values and tilt angles) in A positive epochs. 

(3) The record high cosmic ray intensities observed in 2009 resulted from a reduction in B 
(from ~5.4 nT to —3.9 nT) in comparison with the 1986 minimum, rather than from 
a reduction in the tilt angle (which actually increased slightly in 2009, to ~14°, vs. 
-20° in 1986). 

(4) The tilt angle (a) evolution may be systematically different in even- and odd-numbered 
cycles. In cycles 20 and 22, a appeared to decay more quickly following solar cycle 
maximum than it did in cycles 2 1 and 23. 

(5) The Sun’s message of changing magnetic field is transmitted to the heliosphere pri- 
marily by CMEs at solar maximum and by HSSs at solar minimum. 

(6) The cosmic ray intensity is anti-correlated ( r = —0.82) with the contribution of CMEs 
to solar wind B (based on annual averages), weakly anti-correlated (r = —0.32) with 
the contribution of HSSs to B, and essentially uncorrelated with the SSW contribution. 

(7) The ~30-40% increase in the modulation amplitude of ~2.5 GV cosmic ray electrons 
during the recent cycle minimum suggests that the intensity of ~2 GV protons should 
have increased by a comparable amount (Heber et al. 2009). However, such an increase 
for protons would have been inconsistent with a close relationship between the Newark 
cosmic ray intensity and solar wind B from 1976-2009. Moreover, based on the re- 
lationship between the Newark neutron monitor count rate for this same interval and 
the tilt angle of the HCS, it is difficult to attribute the observed absence of a ~35% 
increase in the ~2 GV proton count rate to a larger tilt angle. 

(8) Long-term reconstructions of solar wind B , both those based on geomagnetic data and 
those inferred from cosmogenic nuclei data from tree rings and ice cores, are in good 
agreement for the 20th century. 

(9) The close relationship between solar wind B and GCR intensity can be used to re- 
produce the modulation curve with reasonable fidelity during the space age, with the 
exception of cycle 20. This relationship can then be combined with the validated long- 
term reconstruction of B to confidently extend the modulation curve back to at least 
M900. 

(10) A recently proposed “floor” in solar wind B of ~2.8 nT implies a minimum (per- 
manent baseline) modulation potential of MOO MV that results in a ceiling in the 
cosmic ray intensity which lies below the local interstellar (LIS) spectrum at any given 
energy/rigidity. The notion of a floor in B is at variance with a recent 10 Be-based re- 
construction (Steinhilber et al. 2010) that shows several episodes when B was ~0 nT 
during the past MO 4 years. 

3.2 Discussion 

While it is generally accepted that barrier modulation is dominant at solar maximum, drift 
effects arc thought to control modulation away from maximum. The solar minimum year 
of 2009 provided a controlled experiment on this point. During 2009, the tilt angle was 
essentially constant in comparison with previous minima while solar wind B dropped, and 
the ~2 GV cosmic ray intensity rose, significantly. Thus it appears that diffusion, driven by 
the minimum-to-minimum variation in B, is the primary modulation process at this phase of 
the solar cycle as well. This picture is supported by the multi-valued relationship between 
GCR intensity and the tilt angle and the more nearly single-valued relationship between 
GCR intensity and B, for both polarity epochs, from 1976-prcsent. 

The Sun sends its message of increased magnetic field (reduced diffusion coefficient) 
to the heliosphere primarily via CMEs at maximum, with high-speed streams playing an 
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increasing role as the cycle declines (e.g., during the mini-cycles in 1974, 1984, and 1994). 
At solar minima, when CMEs occur infrequently, HSSs from polar coronal holes (Wang and 
Sheeley 1994) are the main source of cycle-to-cycle variation in B that drives modulation 
at these times. The high Newark neutron monitor counting rate in 2009, compared to the 
average rate for 1976, 1986, and 1996, resulted from the drop in the average contribution 
of HSSs to B from ~2.8 nT in 1976/1986/1996 to ~1.2 nT in 2009 (thus accounting for 
the drop in B from ~5.4 nT to ~3.9 nT). The drop in 5 hss in 2009 is a direct reflection of 
the ~45% decrease in the solar polar field strength between the minima preceding Cycles 
23 and 24 (Svalgaard et al. 2005; Cliver and Ling 2010). Cliver and Ling (2010) equated 
the SSW with a floor, or constant baseline level, in the solar wind. Such a floor will not 
contribute to either the 1 1 -yr or long-term variation of GCR intensity. 

The effects of large-scale curvature and gradient drifts are most notable at the onset of the 
modulation cycle in A positive epochs (the rise of odd-numbered solar cycles) when GCRs 
respond weakly to increases in B (and tilt angle). The weaker response of the GCR intensity 
to changes in B at such times is attributed to the drift-induced preference for positively 
charged particles to approach the inner heliosphere from the poles at these times (Jokipii 
and Thomas 1981) and the relative confinement of CMEs to low latitudes at the onset of 
the solar (modulation) cycle (Hundhauscn 1993; Gopalswamy et al. 2010). As the cycle 
develops, CMEs form at progressively higher latitudes. Large broad CMEs originating at 
latitudes above ~30-40° can encompass the polar regions and begin to “close the shell” 
in the inner heliosphere (McDonald et al. 1993; Cliver et al. 1993; Cliver and Ling 2001a; 
Lara et al. 2005). As demonstrated by Potgieter and Le Roux (1994), particle inflow from 
over the poles of the Sun during A positive cycles is more rapid than inflow along the wavy 
neutral sheet during the corresponding phase of A negative cycles. Thus GCR recovery times 
are shorter during A positive cycles (Wibberenz et al. 2002). As a result, the onset of strong 
modulation is delayed at the onset of odd-numbered solar cycles until the polar regions are 
affected. 

Thus far during the current cycle 24, the GCR intensity vs. tilt angle relationship in 
Fig. 5(a) deviates significantly from that observed two cycles earlier on the rise of cycle 22. 
The unusually rapid rise of the tilt angle in 2010 is likely related to the weaker polar fields 
of the recent minimum (Svalgaard et al. 2005; Wang et al. 2009), which are eroding as the 
new cycle develops and the equatorial dipole gains strength. Note that despite the recent 
deviation in the GCR intensity vs. a relationship, the relationship between solar wind B and 
GCR intensity (Figs. 3 and 5(b)) is largely unchanged by the data for the onset of the new 
cycle. 

The quadratic relationship in Fig. 10(a) between the GCR intensity and B for both po- 
larity epochs reflects the delayed response of GCR intensity to the onset of modulation in A 
positive cycles. The weaker response of GCRs to changes in low values of B may also be 
due in part to the weaker (relative to CMEs) “cosmic-ray-effectiveness” of the high-speed- 
streams and CIRs (e.g., Burlaga et al. 1991, 1993; Cliver et al. 2003) which dominate the 
variable (non-floor) component of the solar wind at minimum (Cliver and Ling 2010). 

Given the global nature of modulation, it is remarkable that a reconstruction based on 
a single parameter (B) can largely reproduce the observed GCR intensity for the past ~35 
years. The failure to do so during cycle 20 appears to be a relatively rare occurrence. Fig- 
ure 5 in Svalgaard and Cliver (2005) reveals only one other cycle (14) in the last ~ 130 
years where B was similarly disconnected from the sunspot number (and presumably cos- 
mic ray intensity). Our B -based reconstruction of GCR intensity back to 1835 only applies 
to ~2 GV protons. As pointed out by Ferreira and Potgieter (2003, 2004), such a simplified 
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approach does not work at all energies. A compound model involving both drift and diffu- 
sion is required to reproduce modulation curves at lower rigidities, as well as to simulate 
charge-dependent effects at higher energies. 

The recent convergence of reconstructions of solar wind B based on both geomagnetic 
and cosmogenic data gives confidence that we can go further back in time via ice core analy- 
sis. The deep minima (with B ~ 0 nT), which stand out on the long-term reconstruction of 
Steinhilber et al. (2010), raise fundamental questions: (l) Can the solar wind go away? or 
(2) Can an unmagnetized solar wind exist? The revised floor value of ~2.8 nT recently sug- 
gested by Cliver and Ling (2010) accommodates the bulk of the low values in the Steinhilber 
et al. (2010) reconstruction of B for the past ~10 4 years, except for the 9 sharp dropouts in 
which B goes to zero and a comparable number of dips to ~2 nT. A floor in solar wind B 
implies a permanent modulation that results in a ceiling in the cosmic ray intensity below 
the local interstellar spectrum. 

The conundrum of the large increase in the intensity of GCR electrons at the recent 
minimum will need to be resolved. There are three possibilities to explain the discrepancy 
between the electron and proton intensities: (1) the theory is incorrect/incomplete, i.e., elec- 
tron and proton normalized count rates are not equal at minimum; or (2) the observation of 
the electron intensity, or its interpretation (e.g., rigidity value), is in error; or (3) the quanti- 
tative difference in B between the last two minima (~T.2 nT) was large enough to effect a 
qualitative difference in the modulation. 

In closing, we note one point in which all at the meeting were in agreement. The recent 
extension (by subtraction) of the Sun’s dynamic range of activity makes this an exciting time 
to study the Sun and its impact on cosmic rays. 
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